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Very high energy (VHE, energy > 100 GeV) 7-rays undergo pair production with photons 
of the extragalactic background light (EBL). Thus, the intrinsic 7-ray flux of cosmological 
sources is attenuated and the Universe should be opaque to 7-rays above a redshift depen- 
dent energy. Recently, an indication has been found that the Universe is more transparent 
than predicted by a lower-limit EBL model. Here, this indication is confronted with ad- 
ditional VHE 7-ray spectra and different EBL models. Depending on the model for the 
opacity, the indication persist between a ~ 2.6 cr and ~ 4.3 confidence level. 

1 Introduction 

The Universe is opaque for very high energy (VHE, energy > 100 GeV) 7-rays originating from 
cosmological sources above a certain energy. The reason is the interaction of VHE 7-rays with 
low energy photons of the extragalactic background light (EBL), 7vhe 7ebl — > e + e~ [Tj. The 
EBL ranges from optical / ultra-violet (UV) to far infrared wavelengths and has its origin in the 
emission of stars integrated over the history of the Universe as well as star emission absorbed 
and re-emitted by interstellar dust in galaxies [5J for a review] . As a consequence, the intrinsic 
flux of a VHE 7-ray soured^ is dimmed exponentially by exp(— r(z, E)), where r is the so called 
optical depth, an increasing function with the energy E, the redshift z of the source, and the 
EBL photon density. 

Prominent foreground emission makes a direct observation of the EBL difficult [3] , so one 
has to rely on upper limits derived from VHE 7-ray observations [U for recent upper limits], 
lower limits in the optical and infrared inferred from galaxy number counts [5j [6], or EBL 
models 13 [H as examples]. The EBL model of henceforth KD EBL] closely follows the 
lower limits in the infrared and thus predicts a minimal attenuation for VHE 7-rays. It has 
been used in [lOj henceforth HM12] to test if the Universe is even more transparent and indeed 
an indication for a pair production anomaly (PPA) was found. The authors analyze a large 
sample of VHE 7-ray spectra in the following way. All data points of each intrinsic spectrum 
are fitted with analytical functions and the residuals are calculated. For a reasonable fit to the 
spectra, the residuals Xi defined as the diffcrccc between the model flux and the observed flux 
normalized to the statistical uncertainty on the flux (the latter two are absorption corrected) 



1 Most known extragalactic VHE 7-ray sources are blazars, i.e. active galactic nuclei (AGN) with the jet of 
relativistic particles pointed along the line of sight to the observer. 
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should scatter around a mean of zero. If the KD EBL correctly describes the attenuation, than 
this should also be true for the optical thick regime (i.e. r > 2) and this hypothesis is tested 
with a Student's i-test (see HM12 for further details). With this method, HM12 find a ~ 4cr 
indication that the KD EBL model over-corrects the observed data and thus the opacity seems 
to be over-estimated. 

For the present analysis, all spectral listed in Tab. 1 of HM12 are considered together with 
the five extra spectra listed in Tab. [TJ Only the spectrum of 3C 279 has one data point in 
the optical thick regime for the KD EBL. All other spectra only contribute to the (r > 2)- 
distribution if higher EBL density models are considered. Besides the minimal attenuation KD 
EBL, the models of [3 henceforth FRV EBL] and [9] are used to correct the VHE 7-ray spectra 
for absorption. The latter two models predict a very similar attenuation and a low EBL density 
at UV / optical wavelengths resulting in a lower opacity for energies E < 300 GeV than the KD 
EBL. For higher energies, the models result in a stronger attenuation. The effect of these two 
models on the significance of the PPA is, however, non-trivial as it also depends on the actual 
shape of the VHE spectra. It is also allowed for an additional scaling of the optical depth by a 
factor of 1.3. This value is motivated by the detection of the EBL in VHE 7-ray spectra with a 
best-fit scaling factor of ~ 1.3 for the optical depth jTTj . The scaled and unsealed EBL models 
also frame the EBL density allowed by the lower and upper limits. 



Source 


Redshift 


Experiment 


Energy range (TeV) 


Reference 


B3 2247+381 


0.119 


MAGIC 


0.24 - 0.93 


□1 


RXJ0648.7+1516 


0.179 


VERITAS 


0.21 - 0.48 


PI 


RBS 0413 


0.190 


VERITAS 


0.23 - 0.61 


D3 


1ES 0414+009 


0.287 


VERITAS 


0.23 - 0.61 


D22 


3C 279 


0.536 


MAGIC 


0.15 - 0.35 


.16. 



Table 1: List of additional VHE 7-ray spectra added to the analysis conducted in HM12. For 
all other spectra and the corresponding references see HM12, Tab. 1. 



2 Results 

The results of the significances of the i-test for the PPA are shown in Tab. [5J In general, 
the models predicting a higher EBL density result in lower significances. The reason is the 
migration of residuals into the optical thick regime which broadens the distribution and shifts 
it closer to zero (cf. histograms in the top panel in Fig. [T]) However, the overall trend of 
the residual distribution is unchanged (lower panel of Fig. [T]). The black solid line shows a 
smoothed average^ of the residual distribution. For increasing values of r the average also 
increases towards positive values of the residuals. The right most columns of Tab. [5] show the 
results if only spectra are included that have data points in the optical thick regime in KD EBL. 
For all EBL models, these are the spectra with the highest r measurements and all significances 
increase for a unity scaling of r. This underlines that the indication of the PPA is especially 
due to high optical depth measurements. 

The smallest significance, i.e. the least tension between EBL model and data, is found for 
the KD EBL scaled by 1.3. The significance reduces to just above 2 a. However, the high ultra- 
violet / optical density in this scaled model is difficult to reconcile with recent measurements 

2 It is based on a locally weighted polynomial regression (LOESS) 1171 
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of the optical depth with the Fermi-L AT [18]. Using a sample of 150 BL Lac objects up to a 
rcdshift of 1.6, the authors of [18] obtain a best-fit value of the scaling of the optical depth of 
0.90 ± 0.19 for the KD EBL. 







All 


spectra 


Only 


7~kd > 2 spectra 


EBL model 




T x 1.0 


r x 1.3 


T X 1.0 


r x 1.3 


Franccschini ct al. (2008, FRV EBL) 


13 


2.60 


3.50 


3.45 


2.92 


Kneiske & Dole (2010, KD EBL) 


m 


4.34 


2.04 


4.34 


2.61 


Dominguez et al. (2011) 


m 


2.93 


3.51 


3.89 


3.05 



Table 2: Significances (one-sided confidence intervals) of the PPA for the Students t-test for 
different EBL models including either all spectra or only those with data points corresponding 
to t > 2 in the KD EBL (r = r KD )- 



3 Conclusion 

The indication of the pair production anomaly persist with the inclusion of recent VHE 7-ray 
spectra in the sample used in HM12 at a ~ 4.3 a level for the KD EBL. Different EBL models 
that predict a higher EBL density can reduce the significance which is due to the migration of 
data points to the optical thick regime. This shows the natural limitation of the method used 
here: as long as the fit to the spectra is reasonable, choosing a model with a very high EBL 
density will shift most data points into the optical thick regime and the mean of the distribution 
will be close to zero. At the moment, the bulk of VHE emitting blazars has a redshift between 
0.1 and 0.2. Thus, a change in the EBL density at infrared wavelengths has the strongest 
impact on the results. In order to give a definite answer on the significance of the PPA, more 
observations will be necessary, namely of distant objects at several hundreds of GcV, probing 
the ultra-violet / optical regime of the EBL and of sources beyond several tens of TeV to probe 
the far infrared part of the EBL. A similar hint for a low opacity Universe has also been found 
in Fermi-L AT data [15] . 

If the PPA persist, it can be hint for a mechanism that reduces the opacity of the Universe. 
For instance, cosmic rays originating from AGN could also interact with the EBL and produce 
secondary VHE 7-rays that are responsible for the highest energies in VHE spectra [2D]. In 
such scenarios, it is, however, difficult to account for the short term variability in the optical 
thick rcgimeU Another possibility is the oscillations of photons into axion-like particles (ALPs) 
in ambient magnetic fields. For example, it has been shown that the ALP production in intra- 
clustcr magnetic fields followed by a reconversion in the magnetic field of the Milky Way can 
lead to a substantial hardening of VHE spectra [23] . 

Acknowledgments MM would like to thank the state excellence cluster "Connecting Parti- 
cles with the Cosmos" at the university of Hamburg. 

3 The cascade mechanism is especially sensitive to short term variability of distant sources at 2 > 0.2 at 
energies E > 1 TeV which has not been observed so far 1211 . On the other hand, for nearby sources such as 
Markarian421 or Markarian 501, a sizeable amount of cascade emission can only be produced by protons of 
energies E p > 10 19 eV |22| and a short-time variability at energies beyond tens of TeV (i.e. r > 1) would not 
constrain such scenarios. However, another mechanism would be required to explain how 7-ray photons can 
escape pair production in such cases. 
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Figure 1: Distribution of the fit residuals in the different regimes of the optical depth for the 
FRV EBL. Top panel: histograms of the residual distribution for the regimes (from left to 
right) t<1,1<t<2 and 2 < r. The mean, variance, and probability of a Gaussian fit 
are also shown. Bottom panel: residuals as a function of the optical depth. The marker sizes 
correspond to the redshift of the objects and the color coding represents the energies of each 
data point. The black solid line shows a smoothed average and the stars mark the mean of the 
distributions. 
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